Gravitational lensing as a probe of cold dark matter subhalos 
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In the cold dark matter scenario, dark matter halos are assembled hierarchically from smaller 
subunits. Some of these subunits are disrupted during the merging process, whereas others survive 
temporarily in the form of subhalos. A long-standing problem with this picture is that the number 
of subhalos predicted by simulations exceeds the number of luminous dwarf galaxies seen in the the 
vicinity of large galaxies like the Milky Way. Many of the subhalos must therefore have remained 
dark or very faint. If cold dark matter subhalos are as common as predicted, gravitational lensing 
may in principle offer a promising route to detection. In this review, we describe the many ways 
through which lensing by subhalos can manifest itself, and summarize the results from current efforts 
to constrain the properties of cold dark matter subhalos using such effects. 



1. Introduction: Subhalos and satellite galaxies 

The cold dark matter (CDM) model is based on the 
hypothesis that a significant fraction (around 23% [l[) 
of the energy content of the Universe is made up of non- 
baryonic particles, which interact predominantly through 
gravity and have moved with non-relativistic velocities 
since the earliest epochs of structure formation. While 
this scenario has been very successful in explaining the 
formation of large-scale structures in the Universe (galax- 
ies, galaxy groups and galaxy clusters), its predictions on 
subgalactic scales have not yet been observationally con- 
firmed in any convincing way. On the contrary, there are 
at least two features of current CDM simulations that 
appear to be in conflict with empirical data: the exis- 
tence of high-density cusps in the centres of dark matter 
halos, and a rich spectrum of substructures within each 
halo. There is, however, no consensus on how serious 
these problems really are for the CDM paradigm. 

Massive CDM halos are assembled hierarchically from 
smaller halos. As these subunits fall into the potential 
well of larger halos, they suffer tidal stripping of mate- 
rial which ends up in the smooth dark matter compo- 
nent of the halo that swallowed them. Since this is a 
process that may take several billion years to complete, 
many of these smaller halos temporarily survive in the 
form of subhalos (also known as halo substructures or 
subclumps) within the larger halo. According to current 
simulations, around 10% of the virial mass of a Milky 
Way-sized CDM halo should be in the form of subhalos 
at the current epoch [3, y|. Naively, one might expect 
dwarf galaxies to form in these low-mass halos prior to 
merging, which would result in large numbers of satellite 
galaxies within the CDM halo of each large galaxy. A 
long-standing problem with this picture is that the num- 
ber of subhalos predicted by simulations greatly exceeds 
the number of dwarf galaxies seen in the the vicinity_of 
large galaxies like the Milky Way and Andromeda 



This has become known as the "missing satellite prob- 
lem" . A similar lack of dwarf galaxies compared to the 
number of dark halos predicted is also evident within 
group-sized dark matter halos Q. While most studies 
have focused on the discrepancy between the number of 
subhalos and observed satellite galaxies, the discrepancy 
persists in the field population as well. In a sense, the 
missing satellites are just one aspect of a more general 
problem - the mismatch between the low-mass end of 
the dark matter mass function and the luminosity func- 
tion of dwarf galaxies 0] • 

A number of potential solutions to the missing satel- 
lites problem have been suggested in the literature. 
These can be sorted into three different categories, de- 
pending on how they propose to answer the question "Do 
subhalos exist in the numbers predicted by CDM simu- 
lations?" : 

• No. In the first category, we find modifications of 
the properties of dark matter that act to reduce 
the numbers of low-mass halos and subhalos, in- 
cluding warm dark matter Q , self-interacting dark 
matter Q , fuzzy dark matter fioll and dark matter 
in the form of super WIMPs but also models 
of inflation that produce the required cut-off in the 
primordial density fluctuation spectrum 
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observational biases that would put the resulting 
"dark galaxies" 111 [ill o utside the reach of current 
surveys 0, MEl Sj| . While these mecha- 
nisms may be able to solve the missing satellites 
problem as it was originally defined, solutions of 
this type imply that a vast population of low-mass 
CDM subhalos (hosting very faint stellar popula- 
tions or none at all) should still be awaiting discov- 
ery. 

• Yes, but not in our neighbourhood. The final 
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possibility is that the large halo-to-halo scatter in 
subhalo mass fraction may have left the Milky Way 
and Andromeda sitting inside CDM halos with un- 
usually few subhalos compared to the cosmic aver- 



age [2J, |25[ . This would imply that large numbers 



of CDM subhalos (either bright or dark) should be 
awaiting discovery in the vicinity of more distant 
galaxies. 

Gravitational lensing may play an important role in 
the quest to determine which of these different solutions 
is the correct one. If subhalos do exist, lensing can in 
principle be used to detect even those that are too faint to 
be observed through other means. If subhalos do not ex- 
ist, the absence of subhalo-induced lensing effects should 
be able to tell us so. The goal of this review is to explain 
how this can be achieved and to point out some potential 
pitfalls along the way. The material will be described at 
a level comprehensible even for beginning PhD students, 
with focus on the big picture rather than the computa- 
tional details of lensing. 

In section 2, we review the properties of the CDM sub- 
halo population, as inferred by current N-body simula- 
tions. Section 3 outlines the four expected effects of lens- 
ing by CDM subhalos - flux ratio anomalies, astrometric 
effects, small-scale image distortions and time delay ef- 
fects - which are then covered in more detail in sections 
4-7. A number of open questions and future prospects 
are discussed in section 8. 



2. The properties of cold dark matter subhalos 

N-body simulations indicate that the subhalos within 
a galaxy-sized CDM halo follow a mass function of the 
type: 
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with a ~ 1.9 [2, H^l, albeit with non-negligible halo-to- 
halo scatter at the high mass end (M su t, > 5 x 10 s Mq) 
0, Hll . Current simulations of entire galaxy-sized dark 
matter halos can resolve subhalos with masses down to 
M suh - 10 5 M Q , but the mass function may extend 
all the way down to the cut-off in the density fluc- 
tuation power spectrum, which is set by the detailed 
properties of the CDM particles. For many types of 
WIMP s (e.g. neutralinos) this cut-off lies at ~ 10~ 6 A1© 
13, 0, E| ED, HI, but other CDM candidates may al- 
ter this truncation mass considerably. As an example, ax- 
ions may allow the existence of halos with masses as low 
as 1CP 12 Mq 33j), whereas very few halos with masses 
below 10 4 -10 7 Mq are expected in the case of MeV mass 
dark matter [3] • The overall mass contained in resolved 
subhalos (i.e M su b > 10 5 Mq) within a galaxy-sized 
CDM halo amounts to a subhalo mass fraction around 



/sub ~ 0.1, and extrapolating the mass function given 
by eq. JTJ) towards lower masses does not boost this by 
much [2j. 

Since subhalos are more easily disrupted in the central 
regions of their parent halo, the subhalo population tends 
to be less centrally concentrated than the smooth CDM 
component. The spatial distribution of subhalos within 
?"200 ; the radius at which the density of the halo drops 
below 200 times the critical density of the Universe, can 
be described by (35|: 



N(< x) = N(< r 200 ) 
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where x — r/r2oo and N denotes the number of sub- 
halos within a specific radius. It should be noted that 
this result is based on CDM-only simulations, and that 
the presence of baryons within the subhalos may make 
them more resistant to tidal disruption, thereby boosting 
their number densities in the inner regions of their parent 
halos [36] . Some simulations of cluster-mass halos have 
indicated that the spatial distribution of subhalos may 
be a function of subhalo mass, in the sense that high- 
mass subhalos would tend to avoid the central regions 
more than low-mass ones 37, 38|, but this has not been 
confirmed by the latest simulations of galaxy-sized halos 

While the term subhalo is typically used to denote 
clumps located within the virial radius (or, alternatively, 
?"20o) of a large CDM halo, there is also a large number of 
low-mass clumps located just outside this limit [H, |4(| . 
Some of these have previously been bona fide subhalos, 
and others are bound to venture inside the virial radius 
in the near future. Such objects can through projection 
appear close to lines of sight passing through the cen- 
tres of large galaxies, and may therefore be important in 
certain lensing situations. 

The internal structure of subhalos is still a matter of 
much debate. As low-mass halos are accreted by more 
massive ones and become subhalos, substantial mass loss 
occurs, preferentially from their outer regions. The shape 
of the outer part of the subhalo density profile may there- 
fore be seriously affected by stripping, whereas the inner 
profile is left more or less intact. In many lensing studies, 
CDM subhalos are considered to be singular isothermal 
spheres (SIS) or even point masses. This is mainly for 
simplicity - the lensing properties of such objects are 
well-known, but neither observations, theory nor simu- 
lations favour models of this types for the subhalos pre- 
dicted by CDM (see 4l| for references). 

An SIS has a density profile given by: 



Psis 



2nGr 2 ' 



(3) 



where cr v is the line-of-sight velocity dispersion. This 
model has proved to be successful for the massive galaxies 
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responsible for strong lensing (see section 3.1) on arcsec- 
ond scales 0, [43| . The SIS density profile has a steep 
inner slope (p oc r' 3 with (3 = —2), which in the case of 
massive galaxies is believed to be due to the luminous 
baryons residing in their inner regions. This baryonic 
component contributes substantially to the overall mass 
density in the centre, and its formation over cosmological 
time scales may also have caused the CDM halo itself to 
contract, thereby steepening the inner slope of its density 
profile 0, SEE S3- 

Low-mass CDM halos which never formed many stars 
are unlikely to have density profiles this steep. Instead, 
they should resemble the halo density profiles derived 
from CDM-only simulations. The NFW density profile 
EH, with an inner slope of /3 = —1, has for a number 
of years served as the standard density profile for CDM 
halos, and is given by: 



PnfwM = 



(r/r s )(l + r/r s ) 2 



(4) 



where r$ is the characteristic scale radius of the halo and 
Pi is related to the density of the Universe at the time 
of collapse. Modifications of this formula are required 
once halos become subhalos and are tidally stripped. At- 
tempts to quantify the effects of this mass loss on their 
density profile have been made by Hayashi et al. E3 an d 
Kazantzidis et al. [Hot 

Controversy remains, however, over whether the NFW 
profile gives the best representation of CDM halos (and, 
consequently, of subhalos prior to stripping). Based on 
recent high-resolution simulations, some have argued for 
a slightly steeper inner slope (/3 « —1.2; 5l|) with signifi- 
cant halo-to-halo variations, whereas other have favoured 
a far shallower inner slope 0, [52|, 53 , 54j . Inner density 
profiles as steep as that of the SIS model ([3]) are, however, 
unanimously ruled out. While the internal structure of 
subhalos may be relatively unimportant in certain lens- 
ing situations, it can be crucial in others [55|. For lens- 
ing tests that are sensitive to the exact slope of the inner 
density profile of subhalos 4l|, the subhalo-to-subhalo 
scatter in this quantity may also be a very important. 



3. Gravitational lensing by cold dark matter 
subhalos 

The majority of methods aiming to probe CDM sub- 
halos through gravitational lensing focus on detecting 
subhalos outside the Local Group (typically at redshifts 
z « 0.5-1.0). Whereas the subhalos sitting within the 
dark matter halo of the Milky Way will also give rise to 
lensing, these effects would be very difficult to detect and 
separate from other phenomena. It has been suggested 
that CDM subhalos around the Milky Way may be de- 
tectable through the gravitational time delay they impose 
on millisecond pulsars 5(| 57 1, but the time scales and 



probabilities for such events indicate that this is only vi- 
able for the very low-mass end of the subhalo mass func- 
tion (Al < 1 Mq). In what follows, we will therefore 
focus on the lensing situations relevant for subhalos at 
cosmological distances. 



3.1. Strong and weak lensing 

Gravitational lensing effects can be divided into two 
regimes, strong and weak lensing, depending on the level 
of image distortions produced. Strong lensing can occur 
when the surface mass density along a given sightline ex- 
ceeds a certain (redshift-dependent) critical value, and 
is associated with high magnifications, multiple images, 
arcs and rings in the lens plane. In practice, strong lens- 
ing tends to take place when the line of sight from the 
observer to source lies very close to the centre of a galaxy 
or a galaxy cluster, since it is only there that sufficiently 
high surface mass densities are attained. Weak lensing 
occurs in regions of subcritical surface mass density - in 
practice when the lens is located further away from the 
line of sight - and gives rise to small magnifications and 
mild image distortions. The two situations are schemat- 
ically illustrated in Fig. [TJ 

Generally speaking, weak lensing is extremely common 
in the cosmos (at some level, every single light source is 
affected) but inconspicuous, and can only be detected 
statistically by studying a large number of lensed light 
sources. Strong lensing effects, on the other hand, are 
rare but dramatic, and can readily be spotted in individ- 
ual sources. All published strategies for detecting CDM 
subhalos (in the mass range relevant for the missing satel- 
lite problem) through lensing belong to the strong lensing 
category. 



3.2. Macrolensing and millilensing 

Strong lensing can be further divided into subcate- 
gories, depending on the typical angular separation of 
the multiple images produced: macrolensing (> 0.1 arc- 
seconds), millilensing (~ 10~ 3 arcseconds), microlensing, 
(~ 10~ 6 arcseconds), nanolensing (~ 10 -9 arcseconds) 
and so on. When large galaxies (M ~ 10 12 Mq) are 
responsible for the lensing, the image separation typi- 
cally falls in the macrolensing range, whereas individual 
solar-mass stars give image separations in the microlens- 
ing regime. Since all objects with resolved multiple im- 
ages due to gravitational lensing have image separations 
of > 0.1 arcseconds, the term strong lensing is often used 
synonymously with macrolensing. Lensing by objects of 
dwarf-galaxy masses (~ 10 6 -10 10 Mq), like the CDM 
subhalos relevant for the missing satellite problem, have 
been estimated to give rise to millilensing, although the 
exact image separation depends on the internal struc- 
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FIG. 1: Weak and strong lensing. a) Weak lensing occurs when the lens (here illustrated by a gray elliptical galaxy, surrounded 
by a dark matter halo which extends to the outer circle) lies relatively far from the line of sight between the observer (eye) and 
the background light source (star). In this case, where the sightline misses the central galaxy, only a single image is produced, 
subject to mild magnification and distortion. The signatures of this are only detectable in a statistical sense, by studying the 
weak lensing effects on large numbers of background light sources, b) Strong lensing can occur when the dense central region 
of the lens galaxy is well-aligned with the line of sight. The light from the background light source may then reach the observer 
along different paths, corresponding to separate images in the sky. This case is also associated with high magnifications and 
strong image distortions. The angular deflection in this figure, as in all subsequent ones, has been greatly exaggerated for 
clarity. 




FIG. 2: Macrolensed and singly-imaged sources, a) The sightline towards a distant light source passes through many dark 
matter halos with subhalos, but too far from the dense galaxies in the centres of halos for macrolensing to occur. A subhalo 
in one of these halos happens to intersect the line of sight, thereby potentially producing millilensing effects in a singly-imaged 
light source, b) One of the halos happens to lie almost exactly on the line of sight, thereby splitting the background light 
source into separate macroimages. Furthermore, one of the subhalos in the main lens lies on the sightline towards one of the 
macroimages, thereby producing millilensing effects in this macroimage. 



ture of such objects. In this paper, we will for simplicity 
use the term millilensing for all lensing effects associated 
with dwarf-galaxy mass subhalos, regardless of whether 
multiple images are produced or not. The term mesolens- 
ing is sometimes also used to denote this type of lensing 
(with angular scales intermediate between microlensing 
and macrolensing), but since this word also has an al- 



ternative meaning in the gravitational lensing literature 
[He| . we will avoid it here. 
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a) 



b) 





FIG. 3: Small and large sources, a) A galaxy surrounded by a dark matter halo produces multiple images of a small background 
light source (e.g. an optical quasar), b) For a larger background source (e.g. a galaxy or a radio- loud quasar), the macroimages 
may be stretched into arcs or even a complete Einstein ring. 



3.3. Suitable light sources 

In principle, any distant light source may be affected by 
subhalos along the line of sight. The situation is schemat- 
ically illustrated in Fig. [2^. A random line of sight 
towards a light source outside the local volume passes 
within the virial radius of numerous galaxy-sized CDM 
halos (59j . and may hence intersect subhalos anywhere 
along the line of sight. The probability of hitting a sub- 
halo is, however, rather small in this situation, and sight- 
lines of this type may also pass through low-mass field ha- 
los (i.e. the progenitors of subhalos) 6^, 61, H, EU. It 
will therefore be difficult to distinguish between the lens- 
ing effects produced by these two types of lenses. While 
the low-mass end of the field halo population may be very 
interesting in its own right, lensing by such objects is 
more often considered an unwanted "background" when 
attempting to address the missing satellite problem as it 
is currently defined. It may also be hard to distinguish 
complicated intrinsic source structure from extrinsic dis- 
tortions due to lensing (by either subhalos or low-mass 
field halos) in this case. 

Instead, the main targets for attempts to constrain 
the CDM subhalo population using lensing have so far 
been sources that are already known to be macrolensed 
(see Fig. Et>), which in practice means observing either 
multiply-imaged quasars or galaxies lensed into arcs or 
Einstein rings. By doing so, one preselects a sightlinc 
where one knows that there is a massive dark matter 
halo (and supposedly subhalos) located along the line 
of sight. Whether one sees several distinct, point-like 
images, or elongated arcs that approach the form of an 
Einstein ring, mainly depends on the source size: point- 
like sources (quasars in the optical, but potentially also 
supernovae, gamma-ray bursts and their afterglows) give 



distinct images whereas extended sources (galaxies) give 
rise to arcs and rings (see Fig. [3]). The strong magnifica- 
tion produced by the macrolens (large foreground galaxy 
plus dark halo) acts to boost the probability for lensing 
by the subhalo, and typically augments the observable 
consequences of such secondary lensing. In the case of 
multiple images, millilensing by subhalos can moreover 
be distinguished from intrinsic source features, since any 
structure intrinsic to the source should be imprinted in all 
macroimages, whereas millilensing effects will be unique 
to each image. Transient light sources, like supernovae or 
gamma-ray bursts can in principle also be used for this 
endeavour, but no macrolensed sources of this type have 
so far been detected. 

Lensing by subhalos can give rise to a number of ob- 
servable effects, which we describe in the following sec- 
tions: flux ratio anomalies, astrometric effects, small- 
scale structure in macroimages and time-delay effects. 
In what follows, we will focus on subhalos in the mass 
range from globular clusters to dwarf galaxies (~ 10 5 ~ 
10 10 M ), since current predictions indicate that subha- 
los at lower masses may be very difficult to detect through 
lensing effects. 



4. Flux ratio anomalies 

It was noticed quite early that simple, smooth mod- 
els of galaxy lenses usually fit the image positions of 
macrolensed systems well, whereas the magnifications of 
the macroimages are more difficult to explain [ulj . To see 
how this works, a bit of simple lens theory is required. 

Specific relations are expected to apply for the magnifi- 
cations of macroimages close to each other and a critical 
line. Formally, critical lines are the curves in the lens 
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FIG. 4: Different configurations of a four-image lens: a) Fold, b) Cusp and c) Cross. The upper row shows the caustics 
and position of the source (star) in the source plane. The solid line indicates the inner caustic and the dashed line the outer 
caustic. A source positioned inside the inner caustic produces five images. A source positioned between the inner and outer 
caustic produces three images, whereas a source positioned outside the outer caustic will not be multiply imaged. In the case of 
multiple images, one of the images is usually highly de-magnified, so that only four- and two-image lens systems are observed, 
respectively. The lower row shows the corresponding critical lines and resulting observable images in the lens plane. The inner 
caustic maps on the outer critical line and vice versa. A close pair (A, B) and a close triplet (A, B, C) are produced in the fold 
(a) and cusp (b) configurations, respectively. 



plane where the magnification tends to infinity. If crit- 
ical curves are mapped into the source plane, a set of 
caustic curves is obtained. These separate regions in the 
source plane that give rise to different numbers of im- 
ages (see Fig. [4J . The smooth portions of a caustic curve 
are called folds, while the points where two folds meet 
are referred to as cusps. For a background source which 
is close to either a fold (Fig. 2k) or a cusp (Fig- St>) m 
the caustic of a smooth lens, two respectively three close 
images will be produced near the critical line in the lens 
plane. If the source is placed in the center of the caustic, 
the macroimages will form a cross configuration (Fig. 0b). 

All macroimages can furthermore be described as hav- 
ing either positive parity (meaning that the image has 
the same orientation as the source) or negative parity 
(the image is mirror flipped in one dimension relative to 
the source). When taking the image parity into account 
and assigning negative magnifications to negative parity 
images, the sum of the ma gnifi cations of the close im- 



ages should approach zero [65j, l66|, l67|- The following 
relations should then apply for the flux ratio R of a fold 
configuration: 



-Rfold = 



I MA | 



1MB | 



|MA| + Imb I 



0. 



(•5) 



when the separation between the close images (A & B 
in Fig. 2k) is asymptotically small (68J. Here, m repre- 



sents the magnification of a specific image. For the cusp 
configuration (Fig. 2b), the corresponding relation is: 



R c 



Ima| - Imb I + Imc| 



o. 
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1Mb I 



IMC I 



(6) 



However, most observed lensing systems violate these 
relations. This has been interpreted as evidence of small- 
scale structure in the lens on approximately the scale of 
the image separations between the close images. Mag- 
nifications of individual macroimages due to millilensing 
by subhalos would indeed cause the values for -Rfoid & n d 
^cusp to differ from zero fairly independently of the form 
of the rest of the lens [H, M, [zH, [zl 0, H ■ 

A notable problem with this picture is that both 
semi-analytical structure formation models and high- 
resolution ACDM simulations seem to be unable to repro- 
duce the observed flux ratio anomalies, since the surface 
mass density in substructure is lower than that required 



4.1. Complications: Propagation effects and 
microlensing 

Several alternative reasons for the observed flux 
anomalies have been discussed, such as propagation ef- 
fects like absorption, scattering or scintillation in the in- 
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terstellar medium of the lens galaxy 8l| and microlensing 
by stars in the lensing galaxy 82|. Since some sources, 
like quasars, can exhibit intrinsic flux variations on dif- 
ferent timescales, flux ratios may also be difficult to in- 
terpret if the time delay between the macroimages (see 
section 7) is not well known. 

The relevance of propagation effects can be tested by 
supplementary observations of flux ratios at different 
wavelengths, since flux losses due to such mechanisms 
should vary as a function of wavelength. Microlensing 
by stars can be checked for using long-term monitor- 
ing, as this type of lensing is transient and expected to 
introduce extrinsic variability on the order of months. 
Millilensing by halo substructure can on the other hand 
be treated as stationary [7(J. Extended sources (e.g. 
quasars at mid-infrared and radio wavelengths) should 
also be far less affected by microlensing than small, point- 
like sources (quasars in the optical and at X-ray wave- 
lengths). Even though it is often assumed that radio 
observations of quasars are essentially microlensing-free, 
some caution should be applied, since substantial short- 
term microlensing variability is possible in the special 
case of a relativistic radio jet oriented close to the line 
of sight. This phenomenon has been detected in at least 
one multiply- imaged system 



Mid-infrared imaging of lenses is attractive since the 
flux at such wavelengths should be free from differences 
in extinction among the macroimages, in addition to be- 
ing free from microlensing by stars due to the extended 
source size. Such observations can therefore be used to 
test some of the alternative causes for flux ratio anoma- 
lies. Recent studies have used this technique to exam- 
ine several macrolensed quasars with known flux ratio 

anomalies in the optical ia Hi US Hi- 

The mid-infrared 

flux ratios of about half of these systems can be fitted 
with smooth lensing models, which means that only the 
remaining half of these anomalies arc due to millilensing 
by substructures. 



There are also other observational features that argue 
for substructures as the cause of (at least some) flux ra- 
tio anomalies. Negative parity images (so called sad- 
dle images - e.g. the middle image (B) of the close 
triplet in Fig. [4Jd) are often fainter than predicted by 
smooth lens models. This is expected from millilensing, 
as the magnification perturbations induced by substruc- 
ture len sing have been shown to depend on image parity 
ZH 7g, 82, 88]. In contrast, such anomalies cannot be 
attributed to propagation effects since those should sta- 
tistically affect all types of images similarly, regardless of 
their parity. Whether the lensing is due to luminous or 
dark substructures is, however, a different matter. 



4.2. Luminous substructures 

Luminous substructures have been identified in many 
of the lens systems with known flux ratio anomalies. In- 
cluding such substructures in the lens model tends to 
greatly improve the fit to observations. One example of 
such a lens svstem is the radio-loud quadruple quasar 
B2045+265 [89( which exhibits one of the most extreme 
anomalous flux ratios known. Recent deep imaging of 
this system has revealed the presence of a small satellite 
gala xy which is believed to cause the flux ratio anomaly 
[90| ■ Nearly half of the lenses detected in the Cosmic Lens 
All-Sky Survey (CLASS) also display luminous satellite 
gala xies within a few kpc of the primary lensing galaxy 



Recently, there have been studies combining the results 
from simulations and semi-analytical models of galaxy 
formation to investigate if luminous dwarf galaxies might 
be able to explain the frequency of flux ratio anomalies 
observed 9jj, [92J . They find that the fraction of luminous 
satellites in group-sized halos is roughly consistent with 
the observational data within a factor of two, while the 
results for galaxy-sized halos seem too low to explain the 
frequency of luminous satellites within the observed sys- 
tems. The lensing effect of these luminous dwarf galax- 
ies is also somewhat unclear since most satellites found 
in the inner regions of larger galaxies are expected to 
be 'orphan' galaxies stripped of their dark matter halos. 
To investigate this further, higher-resolution simulations 
involving a realistic treatment of the gas processes are 
required. Possible explanations for the discrepancy be- 
tween the expected and observed fraction of luminous 
satellites include dwarf galaxies elsewhere along the line 
of sight mistakenly identified as the lens perturber 
Luminous substructures may moreover be more efficient 
in producing flux ratio anomalies, since they are likely to 
be denser than dark substructures due to baryon cooling 
and condensation [92j |. 

Projection effects are potentially also important for 
flux ratio anomalies caused by completely dark substruc- 
tures as one expects a large amount of line-of-sight struc- 
ture. Although those structures are less effective than 
substructures within the lens galaxy in inducing magni- 
fication perturbations, the overall effect of line-of-sight 
clumps may be significant HE EES. 



5. Astrometric effects 

In macrolensed systems, the presence of halo substruc- 
ture may perturb the angular deflection caused by the 
lens galaxy and thereby the position of macroimages at 
observable levels, so called astrometric perturbations (see 
Fig. E]). 

This method for detecting subhalos has the advan- 
tage of being relatively unaffected by propagation effects 
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FIG. 5: Astrometric perturbations, a) One of the multiple sightlines towards a distant light source passes through a dark 
subhalo. b) The images of the macro-lensed source are observed at the positions of the gray source symbols. Modelling of the 
lens system with a smooth lens potential predicts the position of the upper image at the white source symbol. The subhalo 
close to the sightline of the image causes a deflection on the order of a few tens of milliarcseconds. 



(absorption, scattering or scintillation by the interstellar 
medium) that may contaminate flux ratio measurements. 
Since the astrometric perturbation is a steeper function 
of subhalo mass than flux ratio perturbations, it is mostly 
sensitive to intermediate and high mass substructures 
and therefore probes a distinct part of the subhalo mass 
function [55L Stellar microlensing may complicate 

the interpretation by producing additional shifts of the 
positions of macroimages [951 ] . but such shifts would be 
transient and predominantly affect point-like sources. 

However, the overall size and probability of subhalo- 
induced astrometric perturbations are expected to be 
rather small. Metcalf & Madau [70| used lensing simu- 
lations of random realizations of substructure in regions 
near images and found that it would take substructures 
with masses > 10 8 M that are very closely aligned with 
the images to change the image positions by a few tens of 
milliarcseconds. Such an alignment would be rare in the 
CDM model. Therefore, they suggest to employ lensed 
jets of quasars observed at radio wavelengths, as such 
sources would cover more area on the lens plane. This 
would increase the probability of having a large subhalo 
nearby, but still allow for pronounced distortions due to 
the thinness of the jet. Metcalf [96[ investigated this tech- 
nique further and used it to show that the lens system 
B1152+199 is likely to contain a substructure of mass 
- 10 5 - 1O 7 M . 

Further observational evidence for astrometric per- 
turbations from small scale structure was found in the 
detailed image structures of B2016+112 \$ HI] and 
B0123+437 [99(. In the latter system, a substructure of 
at least ~ 10 6 Af Q would be needed in order to reproduce 
the observed image positions. 

The CDM scenario predicts that there are far more 



low-mass subhalos than high-mass ones (see equation [T]) 
and their summed effect could in principle add up to a 
substantial perturbation. Conversely, since perturbcrs 
positioned on opposite sides around the macrolens gen- 
erate equal but opposite perturbations, the net effect of a 
large number of substructures may cancel out, ensuring 
that rare, massive substructures dominate the position 
perturbation of the images. Chen et al. [94| have in- 
vestigated this by modelling the effects of a wide range 
of subhalo masses and found that all residual distribu- 
tions have very large peak perturbations 10 milliarc- 
seconds). Since the simulation models predict extremely 
few or no substructures in the inner region of the lens, 
the perturbers must be located further away. Therefore, 
it was also inferred that position perturbations of dif- 
ferent images in any lens configuration may be strongly 
correlated. Although these results suggested that rare, 
massive clumps may cause larger perturbations than the 
more abundant smaller clumps, the astrometric pertur- 
bations of the images were considerable even in models 
where no such massive substructures were present. On 
the other hand, these perturbations are at least partly 
degenerate with model parameters of the host halo. 

Since astrometric perturbations are expected to man- 
ifest themselves at (sub-)milliarcsecond levels, high spa- 
tial resolution observations are required which so far are 
mainly achieved by Very Long Baseline Interferometry 
(VLBI) observations of radio-loud quasars. 

However, recent studies have shown that perturbation 
effects of substructure should also be detectable on larger 
scales (~ 0.1 arcseconds) and at shorter wavelengths in 
extended Einstein rings and arcs produced by galaxy- 
galaxy lensing. Peirani et al. [l00[ ] used the perturbative 
method and lens distributions from toy models as well 
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as cosmological simulations to predict the possible sig- 
natures of substructures. They show that when a sub- 
structure is positioned near the critical line, not only as- 
tromctric but also morphological effects, i.e. breaking of 
the image, will occur which are approximately 10 times 
larger and should be easier to detect. 

Other studies have suggested to use non-parametric 
source and lens potential reconstructions to probe small 
perturbations in the lens p otential o f highly magnified 
Einstein rings and arcs (e.g. HqEE^)- Vegetti & Koop- 



mans [103| have used an adaptive-grid method and shown 
that for substructures located on or close to the Einstein 
ring, perturbations with masses > 10 7 M Q respectively 
10 9 Mq can be reconstructed. This technique may then 
be used to constrain the substructure mass fraction and 
their mass-function slope, once a l arger sample of high- 
resolution lenses becomes available [1041 ] . 

With the upcoming generation of telescopes - e.g. the 
Large Synoptic Survey Telescope (LSST), the Joint Dark 
Energy Mission (JDEM), the James Webb Space Tele- 
scope (JWST) and the Atacama Large Millimeter Ar- 
ray (ALMA) - a significant increase in the number a 
macrolensed lenses is expected which will allow the use 
of astrometric perturbations as a probe of the subhalo 
population. 



6. Small-scale structure in macroimages 

When dark objects in the dwarf-galaxy mass range in- 
tersect the line of sight towards distant quasars, image- 
splitting or distortio n on chara cteristic scales of milliarc- 
seconds may occur |l05l Il06j . As already mentioned, 
quasars can currently only be probed on such small scales 
using VLBI techniques at radio wavelengths, but future 
telescopes and instruments may allow similar angular res- 
olution at both optical and X-ray wave lengths 41 1. 



Using VLBI, Wilkinson et al. [1071 ] reported no de- 
tections of millilensing among 300 compact-radio sources 
and was able to impose an upper limit of Q < 0.01 on 
the cosmological density of point-mass objects (i.e., very 
compact objects, like black holes) in the 10 6 -10 8 M© 
range. However, this does not convert into any strong 
limits on the subhalo population, since CDM halos and 
subhalos are not nearly as dense as black holes. Correct- 
ing for this would decrease the expected image separa- 
tions for a millilens of a given mass and yield a probability 
for lensing that is much lower than assumed in their anal- 
ysis. The sources used were moreover not macrolensed 
- this would have made it difficult to make the distinc- 
tion between subhalos and low-mass field halos as the 
main culprits even if any signs of millilensing had been 
detected (see Fig. [2^). 

The effects that a subhalo can have on the internal 
structure of one of the macroimages in a multiply-imaged 
quasar (Fig. are schematically illustrated in Fig. [5] 



For a small, point-like source (e.g. a quasar observed 
at optical wavelengths), the macroimage may split into 
several distinct images with small angular separations 
(Fig. [6H). A larger source (e.g. a quasar at radio wave- 
lengths) may instead exhibit small-scale image distor- 
tions (Fig.[5jD). Even though quasars may display compli- 
cated intrinsic structure when imaged with high spatial 
resolution, such effects can at least in principle be sepa- 
rated from the features imprinted by millilensing, since 
intrinsic structure will be reproduced in all macroimages, 
whereas millilensing effects are unique to each macroim- 
age [127]. 



Yonehara et al. [1081 ] have argued that a significant 
fraction of all macrolensed optical quasars may exhibit 
secondary image-splitting on mil liarcsecon d scales due to 
CDM subhalos. Inoue & Chiba [Toil El^ have explored 
a similar scenario in the case of the extended images ex- 
pected for macrolensed quasars at longer wavelengths, 
and concluded that the small-scale macroimage distor- 
tions produced by CDM subhalos may be detectable with 
upcoming radio facilities such as ALMA or the VLBI 
Space Observatory Programme 2 (VSOP-2). 



The merits of probing CDM subhalos through the 
small-scale structure of macroimages is that, contrary to 
the case for flux ratio anomalies, there is little risk of con- 
fusion due to microlensing by stars or propagation effects 
in the interstellar medium. Glob ular clusters may be able 
to produce similar effects [lllj . and so may luminous 
dwarf galaxies (i.e. the subset of CDM subhalos that 
happen to have experienced substantial star formation), 
but subhalos are expected to outnumber both of these 
populations, at least in most mass intervals. Instead, the 
main problem with this approach seems to be that CDM 
subhalos may not be sufficiently dense to produce mul- 
tiple images on scales that can be resolved by current 
technology. Most studies of these effects have assumed 
that CDM subhalos can be treated as SIS lenses, resulting 
in a gross overprediction of the image separations com- 
pared to more realistic subhalo models 4l| . The angular 
resolution by which macrolensed quasars can be probed 
is on the other hand likely to increase substantially in 
the coming years, in principle reaching w 0.04 milliarc- 
seconds with the VSOP-2 mission (scheduled to launch 
in 2013). 



[127] The distinction between these small-scale changes in the mor- 
phologies of macroimages, and the astrometric effects discussed 
in section 5, becomes somewhat arbitrary in some cases since 
image distortion may both shift the centroid of an image and 
alter its overall appearance (e.g. through the introduction of 
new, small-scale images). The distortion of macrolensed jets is 
for instance usually referred to as an astrometric effect. 
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FIG. 6: A foreground galaxy with a dark matter halo produces multiple macroimages of a background light source. A subhalo 
located in the dark halo intercepts one of these macroimages, which may give rise to a) additional image small-scale splitting 
of the affected macroimage if the source is sufficiently small, or b) a mild distortion in the affected macroimage, if the source 
is large. 



7. Time delay effects 

The images of a macrolensed light source (see Fig. [3^) 
are subject to different time delays, which become de- 
tectable when the source exhibits intrinsic temporal vari- 
ability over observable time scales. These time delays 
stem from a combination of differences in the relativis- 
tic time delays (clocks running slower in deep gravita- 
tional fields, also known as Shapiro time delays) and 
the differences in photon path lengths (due to geomet- 
ric deflection) among the macroimages. Since quasars 
are both non-transient and known to vary significantly in 
brightness on time scales of hours and upwards, they are 
very convenient targets for observing campaigns aiming 
to measure such time delays. At the current time, around 
20 macrolensed quasars have measured time dela ys (w ith 
typical delays of At ~ 0.1-400 days; see Oguri [112j for 
a recent compilation) . Time delays of this type have of- 
ten been used to constrain the Hubble constant and the 
density profile of the macrolens (i.e. the overall gravi- 
tational potential of the lens galaxy and its associated 
dark halo) , but can also potentially be used to probe the 
CDM subhalos of the lens galaxy. 

As shown by Keeton & Moustakas [l!3| . the presence 
of subhalos within the macrolens will perturb the time 
delays predicted by smooth lens models, and may also 
violate the predicted arrival-time ordering of the images. 
Such violations would signal the presence of subhalos in a 
way that, unlike the case for optical flux ratio anomalies, 
cannot be mimicked by dust extinction or microlensing 
by stars. The time delay perturbations due to subha- 
los are typically on the order of a fraction of a day. By 
pushing the uncertainties in the observed time delays to 
this level, strong constraints on CDM subhalo popula- 



tions may potentially be derived. One case of a time 
ordering reversal which may possibly be attributed to 
subhalos has already bee n identifi ed in the macrolensed 
quasar RX J1131 - 1231 [Tl3l [Tl4 . 



If the subhalos themselves give rise to small-scale image 
splittings (as described in section 6), short time lags be- 
tween the light pulses of the separate small-scale images 
would be introduced. This imprints echo-like signatures 
in the overall light curve of astronomical objects with 
short-term variability (such as gamma-ray bursts and X- 
ray quasars), even if the small-scale images cannot be 
spatially resolved. These echos correspond to light sig- 
nals transported through small-scale images with longer 
time delays than the leading image, and the flux ratios 
of the peaks are given by the different magnifications of 
these images. The light curves of gamma-ray bursts have 
been used to search for such light echos in the interval 
~ 1-60 s, resulting in upper limits (fl < 0.1) on p oint- 
mass dark objects in the 10 5 -10 9 M range [l 151 ] and 
even a few can didat e detections of repeating flares due 
to millilensing [lld |. However, just like in the case of 
the searc h for spatial millilensing effects by Wilkinson et 
al. \l0l\ . current investigations of this kind have little 
bearing on CDM subhalos, since the probability for sub- 
halo millilensing is too low when the target objects are 
not macrolensed. Yonehara et al. |l08| instead suggested 
monitoring of macrolensed quasars, predicting that CDM 
subhalos may produce light echos separated by ~ 1000 
s, which could potentially be detected in X-rays, where 
rapid intrinsic flares have been observed. This lensing 
situation is schematically illustrated in Fig. [7] 
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FIG. 7: A galaxy with a dark matter halo produces distinct 
maeroimages of a background light source. If this source dis- 
plays intrinsic variability, observable time delays between the 
different maeroimages may occur. If one of the maeroimages 
experiences further small-scale image splitting due to a sub- 
halo along the line of sight, a light echo may be observable 
in the affected macroimage. This may serve as a signature of 
millilensing in cases where the small-scale images blend into 
one due to insufficient angular resolution of the observations. 



8. Open questions and future prospects 

As we have argued, lensing can in principle be used to 
probe the CDM subhalo population, but has so far not 
resulted in any strong constraints. Most studies have fo- 
cused on flux ratio anomalies, but a number of studies 
now suggest that subhalos by themselves are unable to 
explain this phenomenon [63, 7^, 7^, 8(| 93|. If correct, 
this would limit the usefulness of this diagnostic, since 
some other mechanism must also be affecting the flux ra- 
tios. Luckily, constraints from other techniques, such as 
astrometric perturbations, small-scale image distortions 
and time delay perturbations may be just around the 
corner. 

Observationally, the future for the study of strong 
gravitational lensing is looking bright. As of 2009, around 
200 macrolensed systems have been detected with galax- 
ies acting as the main lens. Planned observational facili- 
ties such as the Square Kilometer Array (SKA) and the 
LOw Frequency ARray for radio astronomy (LOFAR) 
at radio wavelengths and JDEM & LSST in the optical 
have the power to boost this number by orders of magni- 
tude in the coming decade 1171 ]. The spatial resolution 
by which these systems can be studied is also likely to 
become significantly better, approaching ~ 10 milliarc- 
seconds in the optical and ~ 0.1 milliarcseconds at radio 
wavelenghts 41 1. 

On the modelling side, there are still a number of is- 
sues that need to be properly addressed before strong 
constraints on the existence and properties of CDM sub- 



halos can be extracted from such data. 



8.1. Input needed from subhalo simulations 

The largest N-body simulations of galaxy-sized halos 
are now able to resolve CDM subhalos with masses down 
to ~ 10 5 Mq, but there are still a number of aspects of 
the subhalo population that remain poorly quantified and 
could have a significant impact on its lensing signatures: 

• What is the halo-to-halo scatter in the subhalo 
mass function and how does this evolve with red- 
shift? 

• What are the density profiles of subhalos? How 
does this evolve with subhalo mass and subhalo 
position within the parent halo? How large is the 
difference from subhalo to subhalo? 

• What is the spatial distribution of subhalos as a 
function of subhalo mass within the parent halo? 
What is corresponding distribution outside the 
virial radius? 

• How do baryons affect the properties of subha- 
los? Can baryons promote the survival of subhalos 
within the inner regions of their host halos? 

The lensing effects discussed in sections 4-7 are sensi- 
tive to the density profiles and mass function of subhalos, 
albeit to a varying degree 55j . Attempts to quantify the 
effects of different density profile s of lensing signature 
have been made 0, liS El, but the models 
used are still far from realistic, and many of those active 
in this field still cling to SIS profiles for simplicity. 



8.2 The role of other small-scale structure 

CDM subhalos are not the only objects along the line 
of sight to high-redshift light sources that are capable 
of producing millilensing effects. Many large galaxies 
are known to be surrounded by 10 2 -10 3 globular clus- 
ters with masses in the 10 5 -10 6 M Q range. While typ- 
ically less numerous than CDM subhalos in the same 
mass range, they are concentrated within a smaller vol- 
ume (the stellar halo) and have more centrally concen- 
trated density profiles, thereby potentially making them 
more efficient lenses. We also expect a fair share of lumi- 
nous dwarf galaxies within the dark halos of large galax- 
ies. These dwarfs may well represent the subset of CDM 
subhalos inside which baryons were able to collapse and 
form stars, but if so, this means that they may have 
density profiles significantly more centrally concentrated 
than their dark siblings. While the role of globular clus- 
ters and luminous satellite galaxies has been studied in 
the case of flux ratio anomalies 71, 92], their effects on 
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many of the other lensing situations discussed in pre- 
vious sections have not yet been addressed. Low-mass 
halos along the line of sight may also affect these lensing 
signatures, and sometimes appreciably so 62, 63| . 

Aside from dwarf galaxies and globulars, there may 
of course be other surprises hiding in the dark halos of 
galaxies. Intermediate mass (M ~ 10 2 -10 4 Mq) black 
holes, formed either in the very early Universe or as the 
remna nts of populati on III stars may inhabit the halo 
region E3 Ell and could give rise to millilensing 



effects [12l| . If accretion onto such objects is efficient, 
the predicted X-ray properties of such black holes al- 
ready place very strong constraints on their contribution 
to the dark matter (O ;$ 0.005 11221 L but the more gen- 
eralized dynamical 123 . 124 . 125l | and lensing 115l 126| 
constraints on other types of dark objects in the star 
cluster mass range (~ 10 2 -10 5 Mq) are otherwise rather 
weak (fi < 0.1). Lensing observations originally aimed 
to constrain the CDM subhalo population may therefore 
also lead to the detections of completely new types of 
halo substructure. As telescopes attain better sensitivity 
and higher angular resolution in the next decade, we can 
surely look forward to an exciting new era in the study 
of dark matter halos. 
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